Extensivity of Irreversible Current and Stability in Causal Dissipative Hydrodynamics 



G. S. Denicol, T. Kodama, T. Koide and Ph. Mota 

We extended our formulation of causal dissipative hydrodynamics [T. Koide et al, Phys. Rev. 
C75, 034909 (2007)] to be applicable to the ultra-relativistic regime by considering the extensiveness 
of irreversible currents. The new equation has a non-linear term which suppresses the effect of 
viscosity. We found that such a term is necessary to guarantee the positive definiteness of the 
inertia term and stabilize numerical calculations in ultra-relativistic initial conditions. Because of 
the suppression of the viscosity, the behavior of the fluid is more close to that of the ideal fluid. Our 
result is essentially same as that from the extended irreversible thermodynamics, but is different 
from the Israel-Stewart theory. A possible origin of the difference is discussed. 
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It is widely believed that the basic features of collective motion in relativistic heavy-ion collisions can be well 
described by using the (almost) ideal hydrodynamic model Several studies on the effects of viscosity are available 
in various works and seem to support such a vision 0, U 0, H, H, 0, H, S El- However, the inclusion of the 
dissipation into the relativistic hydrodynamics is not trivial because the naive generalization of the Navier-Stokes 
equation gives rise to the problems of acausality and instability. Various theories have been proposed to incorporate 
CI dissipation consistent with causality and stability: the divergence type theory [3| , the Israel-Stewart theory [lB] , the 

! ! extended irreversible thermodynamics [lj|, Carter's theory [l7|, Ottinger-Grmela formulation the approach base 

on conformal field theory [l9[ and so on. However, the relations between these theories are not well clarified and the 
Oh formulation of the relativistic hydrodynamics itself has not yet established. 

One crucial point of these theories is that dissipative fluids behave as a kind of non-Newtonian fluids in a relativistic 
energy regime [l2l . 113 . Recently, we proposed another formulation of the relativistic dissipative hydrodynamics, 
stressing this aspect[8, 9]. The definition of the irreversible current is modified by the introduction of a memory 
l— ~" '■ function, and we showed that such an effect is enough to solve the problem of acausality. We also found that the 
theory is stable against the linear perturbations around the hydrostatic state [l(| • This vision offers new possibilities 
in adapting the techniques developed for Newtonian fluids to a non- Newtonian regime 0, Il2l . [l3| . 

Although our equations of the dissipative hydrodynamics can be reduced to the truncated Israel-Stewart (IS) 
t • . theory, this formulation can be generalized for the application to ultra-relativistic situations. One of the important 
assumptions in the derivation of hydrodynamics is the local equilibrium ansatz: at any space point, there should 
exist an finite extension of the fluid which is described by the thermodynamic laws in equilibrium. We refer such 
an element of the fluid to a fluid cell. The irreversible currents are phenomenologically defined so that the second 
law of thermodynamics is satisfied for each fluid cell. These fluid cells permit superpositions and need not to be 
exclusive, each other, but they have to be of finite size to apply thermodynamics. However, in the usual hydrodynamic 
formulation, the finiteness of the fluid cell is considered to be irrelevant and the irreversible currents are defined by 
applying the second law for densities of extensive quantities. This does not give rise to problems in the (relativistic) 
Navier-Stokes theory, because thermodynamic forces instantaneously produce irreversible currents. That is, the time 
scale of deformation of a fluid cell is considered infinitely large in comparison with this scale. This is also the case 
when the relaxation time, which characterizes the memory effect, is small compared to the variation scale of the fluid. 
However, when the relaxation time is of the same order as the variation scale, we have to distinguish the quantities 
affected by the change in internal degrees of freedom from the global kinematic degrees of freedom in introducing 
memory effects. 

In this paper, we rederive the causal dissipative hydrodynamics by taking the finite size of the fluid cell into 
account. For simplicity, we consider the 1+1 dimensional system. We found that the equation of the bulk viscosity 
has a nonlinear term. Because of the nonlinear term, the effect of bulk viscosity is suppressed and the behavior of 
the fluid is closer to that of the ideal fluid. More importantly, we found that this effect is indispensable to implement 
stable numerical calculations in ultra-relativistic initial conditions. 

This paper is organized as follows. In SecJITl we discuss the new formulation of relativistic dissipative hydrodynamics 
by taking the finite size of the fluid cell into account. The restriction for the parameters which are consistent with 
causality is discussed in Sec lIIII In Sec lIVl we apply our theory to the 1+1 dimensional scaling solution of the Bjorken 
model. In Sec|Vl we implements numerical simulations. We apply the smoothed particle formulation to solve the 
hydrodynamics numerically. Here, we discussed three examples, the shock formation, the expansion to vacuum with 
Landau initial condition and the appearance of nonperiodic oscillations similar to turbulence. We show that the result 
of our formulation can be justified from the viewpoint of the extended irreversible thermodynamics in Sec lVIl The 
relation between our theory and the Israel-Stewart theory is discussed in Sec lVIII Sec lVIIIl is devoted to concluding 
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remarks. 



II. EXTENSIVE MEASURE FOR THE IRREVERSIBLE CURRENT 



For simplicity, we consider the case of vanishing baryon chemical potential for the simple 1+1 dimensional system. 
In this case, the hydrodynamic equations of motion can be written as only the conservation of the energy-momentum 
tensor, 

d^ v = 0, (1) 

together with the thermodynamic relations. We adopt the local equilibrium ansatz in the energy local rest frame, as 
proposed by Landau-Lifshitz |2C| ■ The energy-momentum tensor is given by 

T"" = (e + P + n) u»u" - (P + II) g"", (2) 

where e, P, and IT are, respectively, the energy density, pressure, four-velocity and bulk viscosity. From Eq.(fT]), 
we obtain the entropy production rate in terms of densities, 

3^3^) = ~m^. (3) 

We should remember that the local thermal equilibrium ansatz must be applied to a fluid cell which, in principle, 
has a finite volume defined by the coarse-grained size of microscopic degrees of freedom. The thermodynamic laws 
should be applied to the integrated quantities of the fluid inside each cell. To see this more clearly, let us introduce 
the volume V* of such a fluid cell. That is, V* is the volume of the fluid at the point r, inside of which the fluid is 
considered to be homogeneous and satisfies the thermodynamic laws in equilibrium. The flow of the fluid deforms 
such a cell so that its volume is a function of time. If we follow the fluid flow given by the velocity field v, the time 
variation of V* is given by 

1 dV* 

— = V • v (4 

V* dt w 

or in a covariant form, 

fl„ (tru") = 0, (5) 
where we have introduced the proper reference density a by 

a= v = T*> (6) 

where 7 is the Lorentz factor. 

Now, let us denote the extensive measure of the entropy inside this volume as s — sV = s/a. Then we rewrite 
Eq.© as 

ds 

T— = -JF = -Ud^, (7) 

with 

n 

n = nv = - (8) 

a 

is the extensive measure inside the fluid cell of the irreversible current II. 

One can see that Eq. has the structure that the net entropy production in the cell is given by the product of the 
irreversible displacement II occurred in the cell and the corresponding thermodynamic force field F = d^u^ in the 
cell. It should be noted that, from Eq. ([5]), the thermodynamic force is reexpressed as 



This result means that the thermodynamic force for the fluid cell is given by the change of the cell volume, indicating 
clearly the physical meaning of the bulk viscosity: the resistance to the change of the volume of the system. 
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In the relativistic Navier-Stokes theory (the Landau-Lifshitz theory), it is assumed that the bulk viscosity per 
volume element is produced by the thermodynamic force without any retardation, 

n = -fjF = --a M u". ( 10 ) 

a 

where f] represents a extensive measure of coupling (total charge) for the whole matter inside the cell. We thus 
identify ( as the bulk viscosity coefficient. By multiplying a for both sides, we reproduce the usual result of the 
Landau-Lifshitz theory 20]. That is, the finite size effect does not affect the definition of irreversible currents in the 
Navier-Stokes theory. However, it is by now well known that the derived equation has the problem of acausality and 
instability @, 0, E3l ■ 

To solve these difficulties, we introduce a memory effect to the irreversible current by using a 

memory function. 

When microscopic and macroscopic scales are clearly separated, the time scale of the variation of the fluid cell 
itself is infinitely large compared to the microscopic scales, and only the change in the internal degrees of freedom is 
relevant. If, however, this is not true, they are entangled and we have to distinguish their roles. One is the change 
of the internal degrees of freedom with the transient effects (memory effect) and the other the motion of the global 
kinematic degrees of freedom (fluid cell deformation). As the simplest memory function which can be reduced to the 
differential equation, we apply 

G (r, r>) = — L_ cxp (- f . (11) 



TR (TO r V Jr> ^ 

Here, t r is a relaxation time which characterizes the time scale of the retardation. In the previous work we 
have applied this memory function to the thermodynamic force r/F (not fjF) to introduce the retardation. As a 
consequence, the final form of the derived equation for n is same as that of the so-called truncated form of the IS 
equation, where nonlinear terms are ignored. However as we mentioned before, when we deal with memory effects, 
we should not use densities. This is because memory effects relate different thermodynamic states of the matter for 
different times and this depends on the size of the system not necessarily in an extensive manner. Then the memory 
effect should be applied to the integrated quantity of a fluid cell. 

In this paper, we introduce the memory effect of the thermodynamic force field on an extensive measure for the 
irreversible current. Then the bulk viscosity consistent with causality is given by 

6 (r) = - f dr'G (r, r') ^u" + H exp (- f -^—) , (12) 

J TO <* \ J T„ TROT')/ 

where Ho is the initial value given at tq. This integral form is equivalent to the solution of the following differential 
equation, 



7n^-n + n = -ia M ti' t . (13) 

or a 



This equation can be reexpressed in terms of the density as follows, 



r R ^-n + n = -(c + T R n)d^. (14) 

ar 



On the other hand, the corresponding equation in our previous paper [8(] where the memory effect is directly applied 
to the densities has the form. 



T R -^-n + n = -c^u". (15) 

ar 



The difference of Eq. fTi)) from Eq. (fT5")) is the presence of the term —tbTW^u 1 *. 

It should be noted that Eq. (fT4]) also does not explicitly depend on the cell volume 1/a. This is because the size of 
the cell volume is irrelevant as far as the length is much larger than the mean-free path and much smaller than the 
typical hydrodynamic scale. For later convenience, we call the causal dissipative hydrodynamics without the finite 
size effect, Eq. (fT5")) as "linear causal dissipative hydrodynamics (LCDH)", whereas the one with finite size effect, Eq. 
(fT4|) as "nonlinear causal dissipative hydrodynamics (NLCDH)" . 

An important effect due to the nonlinear term is the lower bound of the bulk viscosity. In LCDH, the bulk viscosity 
can, in principle, take any negative value. In NLCDH, when the bulk viscosity becomes negatively large, the effective 
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bulk viscosity coefficient £ e ff = C + T R.n eventually changes sign and the bulk viscosity start to increase. Thus, the 
bulk viscosity in NLCDH cannot be smaller than 

n min = -— , (16) 

when the initial value of II is larger than n m in - This aspect plays an important role for the stability of numerical 
simulations of ultra-relativistic cases, as we will see Sec. [V] 



III. PROPAGATION SPEED OF SOUND 



We parametrize the bulk viscosity coefficient and the relaxation time as follows, 

C = as, (17) 

TR = FTp 6 ' (18) 

where a and b are arbitrary constants. 

As was pointed out in [§, 0, ll~0j ] . LCDH can be acausal depending on the choice of the parameters. To see the 
limitation of the parameters, we have to calculate the propagation speed of NLCDH. Following [1,0, [13], we discuss 
the linear perturbation around the hydrostatic state. Then the nonlinear term of the equation of the bulk viscosity 
(fT4"]) disappears and the dispersion relation is same as that of LCDH [UGH. Thus, by assuming the group velocity 
gives the propagating speed of the dissipative fluid, we found 

v c = y/l/b + a, (19) 

where a = dP/de. To satisfy causality v c < 1, the parameter b should be less than 1/ (1 — a). This is completely the 
same restriction as the case of LCDH [1, El. 
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FIG. 2: The time evolution of the bulk viscosity of the scaling solution. The dashed and solid lines corresponds to LCDH and 
NLCDH, respectively. 



IV. SCALING SOLUTION 

We apply NLCDH to the one dimensional scaling solution of the Bjorken model. Then the hydrodynamic equations 
are given by 

e + P + n 

d r e + = 0, (20) 

T 

T R a T n + n = -^5 (21) 

T 

where t = rcoshy and x — rsinhy. 

We adopt the massless ideal gas equation of state where a = 1/3. To satisfy causality, the parameter b should 
be larger than 3/2. In this calculation, we choose 6 = 6. For the initial condition, we set e (tq) = 1 GeV/fm and 
n (to) = at the initial proper time tq = 0.1 fm. 

In Fig. [U we show the energy density e as a function of the proper time r. The dotted, dashed and solid lines 
correspond to the ideal hydrodynamics, LCDH and NLCDH, respectively. Because of the memory effect, the behaviors 
of LCDH and NLCDH are similar to that of the ideal fluid at the early stage of the time evolution. After the time 
larger than the relaxation time, the behaviors of LCDH deviates from that of the ideal hydrodynamics. On the other 
hand, the behavior of NLCDH stays close to that of the ideal hydrodynamics, that is, the effect of the bulk viscosity 
is suppressed in NLCDH compared to LCDH. This is directly observed from the behavior of the bulk viscosity as is 
shown in Fig. [2] 

It should be noted that the scaling solution of the bulk viscosity including the similar nonlinear term was already 
discussed in Q, where the equation derived by Israel and Stewart is discussed. The relation of our equation and the 
Israel-Stewart equation will be discussed in Sec. IVIII 



V. NUMERICAL SIMULATIONS 
A. Smoothed particle formulation 



To solve numerically the hydrodynamic equations, we use the Smoothed Particle Hydrodynamic (SPH) method. 
The original idea of the SPH method is to obtain an approximate solution of hydrodynamics by parameterizing the 



fluid into a set of effective particles [22j. Because of its flexibility to adapt to complex geometries, the SPH method 



has also been extensively applied to the relativistic heavy ion reactions to perfom an event-by-event analysis of the 
data [2l|. For the sake of convenience, we reproduce below the ideas and basic equations of the SPH approach shown 
in @. 
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Let us consider a distribution a (r, t) of any extensive physical quantity, A. In a system like the hot and dense 
matter created in heavy ion collisions, the behavior of a (r, t) contains the effects of whole microscopic degrees of 
freedom. We are not interested in the extremely short wavelength behavior of a (r, t) but rather in global behaviors 
which are related directly to the experimental observables. Therefore, we would like to introduce a coarse-graining 
procedure for a. To do this, we introduce the kernel function W (r — f , h) which maps the original distribution a to 
a coarse-grained version ace as, 

<z C G(r,i) = J a(r,t)W(r-r,h)dr (22) 



where W is normalized, 



J W (f , h) dr = 1, (23) 



and has a bounded support of the scale of h, 

W{r,h)^0, \r\>h, (24) 

satisfying 

IimW(r,fc) =S(r). (25) 

Here, h is a typical length scale for the coarse-graining in the sense that the kernel function W introduces a cut-off in 
short wavelength of the order of h. Thus we will take this value as the scale of coarse graining in the QCD dynamics 
(i.e., the mean- free path of partons) to obtain the hydrodynamics of QGP (h ~ 0.1 fm). 

The second step is to approximate this coarse-grained distribution ace ( r ; t) by replacing the integral in Eq. (|2"2")l by 
a summation over a finite and discrete set of points, {r a (t), a — 1, ... , -/Vsph}, 

JVsPH 

a S P H (r,i)= A a (t)W(\r-r a (t)\). (26) 

a=l 

If the choice of {A a (t), a = 1, . . . ,JVsph} and {r a (t), a = 1, . . . , Nsph} are appropriate, the above expression 
should converge to the coarse-grained distribution aco for large ./Vsph- Parameters {A a (t),a = 1, . . . , A^sph} and 
{r a (t),a = 1, . . . , A*sph} should be determined from the dynamics of the system. In practice, we first choose the 
reference density a* which is conserved, 

^ + V.j = 0, (27) 
where j is the current associated with the density a* . Then, we note that the following ansatzes, 

^sphM)= ]T u a W(\T-r a (t)\), (28) 

a=l 

Ns-ph , /,\ 

j S P H (r,t)= ]T u a ^lW(\r-r a (t)\), (29) 

a — 1 

satisfies the equation 



dt 

where 2/ a ' s are constant. By using the normalization of W, we have 

-'Vsph 



^ PH -V-j S P H = 0, (30) 



/ a*(r,t)d 3 r = V v a . (31) 
-'sph a=1 

Then we can interpret the quantity v a as the conserved quantity attached at the point r = r a (t). Therefore, the 
distribution erg PH (r, t) is a sum of small piece-wise distribution, carrying the density, 

v a W{\v-v a {t)\). (32) 
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These pieces are referred to as " SPH-particles" . 

Using the above reference density and the extensive nature of A, we can write A a in Eq. tpfi]) as 

A a {t) = v a a (33) 

which represents the quantity A carried by the SPH particle at the position r — r Q (t). In fact, the total amount of A 
of the system at the instant t is given by 

Nsph 

A(t)=J2 A «( t )- ( 34 ) 

a=l 

In the ideal fluid, the entropy density is chosen as the reference density and the dynamics of the parameters 
{r a (t), a = 1, . . . , -/Vgpn} are determined from the variational principle from the action of ideal hydrodynamics. The 
entropy density is, however, not conserved for the dissipative fluid. Thus we introduce a new conserved quantity, the 
specific proper density a, which is defined by the flow of the fluid, 

(on") = 0, (35) 

and we will use it as the reference density for viscous fluids. Here, the four-velocity it" is defined in terms of the local 
rest frame of the energy flow (Landau frame) . The specific density is expressed in the SPH form as 

Nsph 

a*(T,t)=Y t v a W(\T-T a {t)\), (36) 

a=l 

where a* = era is the specific density in the laboratory frame and v a is the inverse of the specific volume of the SPH 
particle a. In this work, the specific volume should be interpreted as the volume of the fluid cell, and hence v a is the 
inverse of the cell volume. However, as we showed, the finial results do not depend on this choice and we set v a = 1 
for simplicity. As for the kernel W (r), we use the spline function. 

It should be mentioned that this procedure is only possible provided that the lines of flow in space defined by 
the velocity field it" do not cross each other during the evolution in time. That is, if there appear turbulence or 
singularities in the flow lines, the above definition of Lagrange coordinates can fail. 

Now we apply this method to NLCDH in 1 + 1 dimension. We have to solve the evolution equation of the bulk 
viscosity in the SPH scheme. For this, we express the viscosity as 



n= V v a - M/(|r-r Q (t)|), (37) 




The time evolution of the term II Q can be calculated as 

^T = -- r-(*W) a -— ft, (38) 

where j a is the Lorentz factor of the a-th particle. In the following, we denote the quantity in the observable frame 
with the asterisk. It should be noted that we solve Eq. (fT3"|) instead of (fTi|) in the numerical calculations. 
At the same time, using the SPH expression for the entropy density s* in the observable frame, 



WsPH 

E 

a=l 



v a (-) W{\v-v a {t)\), (39) 



and the evolution of the entropy per SPH particle is given by 

!(;).--?!<*■".• 

where s — s* /u is the proper entropy density. In the above expressions, the relaxation time tr, viscosity coefficient £ 
and temperature T are functions of space and time, so that they should be evaluated at the position of each particle 
o. 
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Finally, we need to express the momentum conservation equation by the SPH variables. We write the space 
component of energy-momentum equation of continuity in terms of the reference density, 



dr 



e + P + U 



d, (p + n) = o. 



(41) 



It should be noted that there exist ambiguities within the resolution of the coarse-graining size h to express the 
equation of motion in the SPH form. However, in the ideal fluid, the SPH equation of motion can be derived by 
the variational method uniquely. Thus, we obtain the equation of motion by using the same SPH parametrization to 
Eq.@IJ|, 



P n + IL 



dr a 



0=1 



( 



Pp + Upap P a +U 



V 



diW(\r a -Tp(t)\) 



(42) 



where the right hand side of Eq.((42|) corresponds to the term di (P + U) written in terms of the SPH parametrization. 
We remark that in the case of vanishing viscosity our result is reduced to the expression derived with variational 
principle for ideal fluids. 

By separating the acceleration and force terms in Eq. (|42[) , we obtain our final expression of the equation of motion 
for each SPH particle, 



where the mass matrix M and the force term F are defined as 



with 



Mij = 7 (c + P + n) Si. 
F j =-d j (P + Il)+Bu j , 

C 



Au 



A = - 
B = A 



1 



7 

>f_ da_ 
a*~dl 



a(e + P + n) + — + n 



n 

tr' 



where a = dP/de. It should be noted that the expression of A in NLCDH is different from that of in LCDH @, 



^lcdh = — - 



a (e + P + n) 



TR 



(43) 



(44) 
(45) 



(46) 
(47) 

(48) 



To carry out the calculation in this scheme, the mass matrix M must be nonsingular. However, this is not guaranteed 
in LCDH. For example, consider the 1+1 dimensional system with ultra-relativistic fluid velocity, u « 7. Then the 
mass in LCDH becomes 



Mlcdh ~ 7 



(l-a)(e + P + II) 



TR 



(49) 



As we discussed the bulk viscosity in LCDH does not have a lower bound, we see that Mlcdh can be zero, while in 
NLCDH, we get 



Mnlc 



DH 



7 



(l-a)(e + P + II) 
7(1 -a) (e + P) [1 



n 




I) 




> 0. 







>7(l-a)(e + P- — 
tr 



(50) 



Here, we used the expression of the bulk viscosity coefficient and the relaxation time, Eqs. |TTJ) and (Tj"5|) . The mass 
does not vanish in NLCDH and the simulation is stable even for the ultra-relativistic situations. 
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FIG. 3: The temperature in the shock formation calcu- FIG - 4: The velocity m the shock formation calculated 
lated in LCDH with a = 1 at t = 0.75 fm, starting from m LCDH wlth a = 1 at t = °- 75 fm - The mltlal velocity 
the homogeneous initial condition (dotted line). ( dotted llne ) at the maximum is 7 = 5. 





FIG. 5: The mass matrix in the shock formation calcu- 
lated in LCDH with a — 1 at t = 0.75 fm. The mass 
matrix crosses zero around the minimum. 



B. Shock formation 



As was discussed in Q, in the LCDH scheme, the numerical calculation can be carried out with the help of 
the additional viscosity for 7 = 2. Here, we show that this scheme becomes unstable for ultra-relativistic initial 
conditions. We use the same additional viscosity as the one proposed in [{|. In Figs. [3] and [4] we show, respectively, 
the temperature and velocity profiles calculated in LCDH with a = 1 at t = 0.75 fm for the initial triger velocity 
7 = 5. The dotted lines are for the initial condition. We can see that the calculation becomes unstable and rapid 
oscillations appear around x — and the calculation eventually collapses. This is due to the vanishing of the mass 
term. The behavior of the mass term is shown in Fig. [5l We can see that the mass becomes zero at the point where 
the calculation shows the rapid oscillation. 

On the other hand, as was shown in the previous section, the mass does not vanish in NLCDH. Thus we can 
implement the numerical calculation with the ultra-relativistic initial condition without numerical singularities. In 
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FIG. 6: The temperature in the shock formation calcu- 
lated in NLCDH with a = 1 and b = 6 at t = 0.75 fm, 
starting from the homogeneous initial condition (dotted 
line) . 



FIG. 7: The velocity in the shock formation calculated in 
NLCDH with a = 1 and b = 6 at t = 0.75 fm. The initial 
velocity (dotted line) at the maximum is 7 = 5. 




FIG. 8: The mass matrix in the shock formation calcu- 
lated in NLCDH with a = 1 and b = 6 at t = 0.75 fm. 
The mass matrix does not cross zero. 



Figs. [S] and El we show, respectively, the temperature and velocity distributions calculated in NLCDH with a = 1 
and b = 6 at t = 0.75 fm, starting from the same initial condition. We can see that the calculation remains stable, 
because the mass does not vanish as is shown in FigjS] 

In all simulations of the present work, we consistently use the additional viscosity in the NLCDH scheme. That is, 
we introduce the nonlinear term for the coarse-grain viscosity used in Q . 



C. Landau initial condition 



Here, we discuss the expansion of the fluid to vacuum and compare the cooling process of LCDH and NLCDH. We 
use the Landau initial condition where the initial temperature is 590 MeV and the initial size is 0.7 fm. In Fig[9j we 
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FIG. 9: The evolution of the temperature with the Lan- 
dau initial condition using a = 0.1 for t = 1, 2 and 3 
fm. The solid and dotted lines denote the calculations in 
NLCDH and LCDH, respectively. 
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FIG. 10: The evolution of the bulk viscosity with the 
Landau initial condition using a — 0.1 for t = 1, 2 and 3 
fm. The solid and dotted lines denote the calculations in 
NLCDH and LCDH, respectively. 




FIG. 11: The evolution of the pressure (dotted line) and 
bulk viscosity (solid line) in NLCDH using a — 0.1 for 
t = 1, 2 and 3 fm. 



show the evolution of the temperature with a — 0.1 and 6 = 6 for t = 1, 2 and 4 fm, from the top. The solid and 
dotted lines represents the results of NLCDH and LCDH, respectively. One can see that the cooling and expansion 
of the fluid of the NLCDH are faster than that of LCDH, similarly to the case of the scaling solution. This is because 
the bulk viscosity in NLCDH is supressed by the nonlinear term in comparison with LCDH. This is explicitly shown 
in Fig. rjUl where the evolutions of the bulk viscosity are plotted. 

As was discussed in Q, the propagation to vacuum in LCDH, a stationary wave is formed and the pressure and 
the bulk viscosity should satisfy the relation P = — H at the boundary. In Fig. [TlJ the pressure (dotted line) and 
the bulk viscosity (solid line) of NLCDH are shown. One can see that the relation P = — H at the boundary is still 
satisfied even in NLCDH. 
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FIG. 12: The evolution of the temperature with the Lan- 
dau initial condition in NLCDH using a = 1 and 6 = 6 
for t = 1.44, 2.04 and 2.64 fm from the top. 



FIG. 13: The trajectory of the fluid element at the central 
rapidity region on the phase diagram as a function of f 
and Ro for a = 0.5, 0.8, 1 2, 3, 4 and 5 from the top fixing 
6 = 6. 



D. Nonperiodic oscillations 

In [l(| , we pointed out that the numerical calculation of LCDH becomes unstable near the central rapidity region 
and the nonperiodic oscillations appear. It is interesting to note that this scenario persists even in the NLCDH 
scheme. In Fig|T2l we show the evolution of the temperature calculated in NLCDH with a = 1, b = 6 for t = 1.44, 
2.44 and 2.64 fm from the top. We use the Landau initial condition with the initial temperature 590 MeV and the 
initial size 0.7 fm. One can see that the nonperiodic oscillations evolve with time in the center of the fluid. 

Interestingly, the appearance of nonperiodic oscillations has a regularity. We investigate the parameter dependence 
of how these oscillations emerge. In FigfT3l we plot the trajectories of the fluid element at the central rapidity region 
as a function of f and R, for the same initial condition. Here, f is the scaled proper time f = t/tr and R is the 
Reynolds number defined as R = — (e + P) /H [l(|. The dotted lines show the trajectories for a — 0.5, 0.8, 1, 2, 3, 4 
and 5 from the top, respectively. The parameter b is fixed to 6. We follow the trajectory of the fluid element in this 
plane up to the point where the oscillation emerges for each value of a. Thus, the line formed by the endpoints of 
these trajectories defines the critical line for the appearance of the oscillations, indicated by the solid line. One note 
that the trajectories without oscillation, a — 0.5 and 0.8, do not cross this line. We confirm numerically that, for 
various sets of parameters, the oscillations appear only when a trajectory crosses the critical line. Instabilities have 
also been analised for the scaling solution in the framework of the first order theory [TlT ] . 



VI. EXTENDED IRREVERSIBLE THERMODYNAMICS 



In this paper, we derived our equation by introducing the memory effect and the finite size effect. As was pointed 
out, there are several different approaches to derive the relativistic hydrodynamics consistent with causality. In this 
section, we briefly review the derivation based on the extended irreversible thermodynamics (EIT) [16| and show that 
our formulation and the extended irreversible thermodynamics give same conclusion. 

The usual thermodynamics describes the thermal equilibrium state which can be described by the so-called ther- 
modynamic variables; energy, volume and number of particles. The extended irreversible thermodynamics is the 
extension of the usual thermodynamics so as to describe the non-equilibrium state, which is characterized by not only 
the thermodynamic variables but also irreversible currents. Then the first law of thermodynamics in the local rest 
frame is 



(51) 
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where EC, 7r y and v % are, respectively, the bulk viscosity, the shear viscosity and the heat conduction times V = l/cr, 
following the definition of this paper. On the other hand, the entropy is expanded around the equilibrium state of 
(E, V, N) as 



s(E,V,N,U,Tr,q) -S(E,V,N\ 



1 d 2 S 



2 SIP 



IT + 



i d 2 s 



2 dTT l idfr lm 



7T J 7T 



1 d 2 s 



2 dv l dvi 



(52) 



where the suffix denotes 



0. It should be noted that we omitted the mixed derivatives with different 



tensority like d 2 S / dlidv 1 , for simplicity. As a matter of fact, such a term can be small according to the Curie principle 
[23| . Here, we used that the first derivatives vanishes because the entropy should be maximum in the equilibrium, 
and ignored higher order derivatives. By comparing the two expressions, we have 



dS n 1 



dS 
dS_ 



i 

'f 



(53) 
(54) 
(55) 



where fii/T is a function only of the usual thermodynamic variables, (E, V, N). Inserting these results into the first 
law (ISTl). we obtain 



TdS = dE + PdV - [idN ~ /3iIMI - /3 2 tt ^dn^ - fcv„dv», 



(56) 



expressed in a covariant form. From the equation of continuity of the energy-momentum tensor and particle flux, we 
have 



a— e + (P + II) + u,d v Ti^ = 0, 
dr 

a^-h + d u ^ = 0. 
dr 



(57) 
(58) 



By combining Eqs. ([56]) . (|57|) and (|58[) together, we have 



(59) 



where 



n 



T dr ' M 

and the entropy four flux is defined by 



T dr T I T 



d^ v 



dr 



(60) 



T 



(61) 



To satisfy the algebraic positivity of this entropy production, we obtain the equations of the irreversible current, 



d^ 1 + Piv^r- = -aiTL, 



dr 



P» vap [d a u fs ~f3 2 a 



dr 



T T dr 



(62) 
(63) 
(64) 



where the projection operators are defined by 



pp,v 



pfiisafi 1_ fpfia pvj3 _|_ pf-if3 pva\ 1 pi±v paf3 



(65) 
(66) 
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with D is the spatial dimension. Here, ctj is a positive parameter. One can easily see that the equation of, for example, 
the bulk viscosity is nothing but the equation obtained in our formulation, by setting £ = 1/ai and tr = f3\ja\. 

To our best knowledge, this is the derivation of the relativistic hydrodynamics based on the extended irreversible 
thermodynamics for the first time. As just described, our formulation and the extended irreversible thermodynamic 
derives the same equation. The only difference is that the thermodynamic variables e and P in NLCDH satisfy usual 
thermodynamic relation, while the variables obey the extended thermodynamic relation in the hydrodynamic equation 
of the extended irreversible thermodynamics. 

Note that it is sometimes said that the concept of thermodynamics is extended in the IS theory. However, as is 
discussed in Appendix in detail, the first law of thermodynamics is modified in a different way in the IS theory, 
without increasing the number of thermodynamical variables. 



As was shown in the previous section, the equation obtained in our formulation can be derived also from the 
extended irreversible thermodynamics. 

Exactly speaking, the similar non-linear term appears even in the IS theory, although in the so-called truncated 
version of IS theory the non-linear terms are ignored. Then the equation of the bulk viscosity in the original IS theory 
is given by 



The last two terms on the r. h. s. are ignored in the truncated IS theory. For detailed derivation, see Appendix [A] 
Here we used the relation £ = l/a± and tr = 0o/ai. 

One can easily see that there are two differences between NLCDH and IS theory. One is the coefficient of the 
non-linear term r^Ild^u^. In NLCDH, the coefficient is given by just 1, but it is 1/2 in the IS theory. The other is 
the last term of Eq. ([57]) . which does not appear in NLCDH. 

As was shown in Sec. IV Al the hydrodynamic equation is stabilized by the non-linear term in NLCDH. In the IS 
theory, we could not show the stability of the IS theory analytically because of the last term of Eq. (|67j) . However, the 
numerical simulation shows that the IS theory is more stable than the truncated IS theory (or LCDH). We calculate 
the shock formation and the fluid expansion with the same parameters and initial conditions, as was discussed in Sec. 
IV Bl and IV CI in the IS theory. We found that the numerical calculations are stable as in the case of NLCDH and the 
behaviors of the IS theory is similar to that of the NLCDH. In Fig[l4l we plot the entropy production as a function of 
time in the calculation of the expansion to vacuum with the Landau initial condition. One can see that the entropy 
production of NLCDH and the IS theory is smaller than that of LCDH (or the truncated IS theory) and the IS theory 
is most close to the Ideal fluid. 

The similar reduction due to non-linear terms also appears even for the shear viscosity as was numerically studied 
in 0]. In the bulk case, we can show explicitly that the suppression happens because of the minimum value of n 
guantced by the nonlinear term. In the case of the full IS theory, it is not obvious why this happens, but surprizingly 
the net effect is very close to ours. 

The non-periodic oscillation appears even in the IS theory. It should, however, be noted that in these simulations, 
we introduced the additional viscosity @. If we do not use the additional viscosity, we cannot implement stable 
numerical simulations even in the IS theory. 

It is important to note that the IS theory also can be derived in our formulation. So far, we employed the 
memory effect between n and F for simplicity reasons. Suppose we apply the same procedure to a/ tr/(T<j H and 
£y/ tr/£T<t d^u^ by using the same memory function, 



VII. RELATION BETWEEN NLCDH AND THE ISRAEL-STEWART THEORY 




(67) 




(68) 



This leads to the following equation for n, 



(69) 



This equation is nothing but the equation of the bulk viscosity in the IS theory (See Appendix [A"| . It can also be 
derived in the framework of the internal- variable theory, although the concept of thermodynamics should be extended. 
See Appendix [5] and [Bl respectively. 
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FIG. 14: The evolution of the entropy production with the Landau initial condition using a = 0.1. The dotted, solid and 
dashed lines represent the calculations of LCDH, NLCDH and the IS theory, respectively. 



VIII. CONCLUDING REMARKS 



macroscopic 


l/V » TR 


> TR 


1/G> M ~ TR 


scale 


(non-relativistic) 


(relativistic) 


(ultra- relativistic) 


irreversible 


linear response 


+ memory effect 


+ finite size effect 


currents 


( J = VF ) 


(J = JG V F) 


(J = fGfjF) 



In this paper, we extended the previous derivation of the causal dissipative hydrodynamics to take into account the 
finiteness of fluid cells and the memory effects on the extensive measure of irreversible current inside the fluid cell. 
The new equation has a non-linear term which suppresses the effect of viscosity. Thus the behavior of this fluid is 
more close to that of the fluid described by the so-called truncated Israel-Stewart theory, where there is no finite size 
effect. More importantly, we found that the non- linear term is necessary to implement stable numerical simulations 
for the ultra-relativistic situations, like high initial velocity, high initial energy density and so on. 

In this study, we found that there are three stages in the structure of hydrodynamics as is summarized in the above 
table. When the time scale of microscopic degrees of freedom are clearly separated from those of the hydrodynamic 
variables, we can assume that the irreversible currents J are immediately produced by the thermodynamic forces F, 
that is, J — rjF. This is realized in the non-relativistic cases, l/d^u^ 3> tr, because the time scale of constituent 
particle of the fluid is much faster than the velocity of the fluid. However, in the relativistic fluids where l/d^u^ > tr, 
the clear separation of the time scales is not necessarily borne out and we have to take into account the retardation 
effect in the formation of the irreversible currents by introducing memory functions G, that is, J = J G rjF. In 
the ultra-relativistic limit where the relaxation time is same order as the scale of inhomogeneity, tr <~ l/d^u 1 *, it is 
important to consider the effect of the finite fluid cell volume 1/cr, because the volume of the fluid cell changes during 
the hysteresis, J = J G fj F. 

We further showed that our formulation and the extended irreversible thermodynamics lead to the same hydrody- 
namic equation. The only difference is that the thermodynamic variables e and P in NLCDH satisfy usual thermo- 
dynamic relation, while the variables follows the extended thermodynamic relation in the hydrodynamic equation in 
the extended irreversible thermodynamics. This may indicate us the robustness of our equation. 
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method 



derived equation 



IVT (see Appendix EJ) 
simplest moment method 
memory function method (not extensive) 



a- 



IS equation 



EIT 

memory function method (extensive) 
equation of this paper 



On the other hand, our result is different from that of the IS theory. We showed that the IS theory also seems 
to be applicable to the ultra-relativistic cases, although we could not show the positivity of the mass explicitly. 
The quantitative difference between our theory and the IS theory is very small but the IS theory is more close to 
the behavior to the ideal fluid. It is also interesting to mention that the equation of the IS theory can be derived 
even in our formulation if we introduce a very peculiar form of thermodynamical variable to introduce the memory 
effect. Although our choice seems more natural and simple from the point of view of the memory function on 
thermodynamical forces, we need experimental and other theoretical supports to decide the appropriate forms of 
thermodynamical variables. The way of extending the thermodynamics to irreversible domain is not unique (See 
also the discussion in Appendix iBl and Ref. [l6|). The schematic mapping of different theories are summarized in 
Fig. IVIIII It should be noted that to derive the IS equation in the memory function method, we have to break the 
extensivity of currents because we employ the memory effect to the quantity proportional to l/\/a as is shown in Eq. 



We are interested in the dynamics of fluid, which cannot be described by the simple Boltzmann equation. It is not 
obvious but if the theory is still applicable to the dynamics of a dilute gas as in the non-relativistic case, the equation 
should be justified from the kinetic argument such as the moment method. As a matter of fact, the problem of the 
simplest moment method is known and there are several proposals for the improvement 25]. The derivation of our 
equaiton from the kinetic point of view is still an open problem. 

T. Koide acknowledges useful discussions with D. Jou. This work is supported by FAPERJ and CNPq. 



APPENDIX A: ISRAEL-STEWART THEORY 



Similar to the extended irreversible thermodynamics, the equations of irreversible currents are derived by applying 
the algebraic positivity of the entropy production in the IS theory. In the derivation of Israel and Stewart, they used 
two equations; one is the definition of the entropy four flux, 



°IS 



s (T, M )^-^-Q", 



and the other is the expression of Q M ? which is assumed by the general quadratic form, 

1 



y (/3 n 2 + 027^7^ + fcv^) 



Here s (T, /i) denotes the entropy density in equilibrium for given temperature T and chemical potential /i. 
By using the algebraic positivity of the entropy production d^S^g > as usual, we obtain 



dn 



2 dr 



a 2 ^ = d, 



0: 



dn, 



—Tr a pd\u - (3 2 , 
2 dr 



a/3 



T d 



a r>uu ( a f* & A a & dVv 1 d ( 



(Al) 

(A2) 

(A3) 
(A4) 
(A5) 



The derivation of Israel and Stewart depends on the validity of the assumption (|A1[) . To derive this equation, it is 
uscfull to derive the following relation in the equilibrium, 



d^ - 



a rpilV _ P N n 
Pv (0) T (0)' 



(A6) 
(A7) 
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here, the subscript (0) indicates quantities in equilibrium and /3„ = u v /T. 

The fundamental assumption used by Israel and Stewart is the so called "release of variations" , which assumes that 
(|A6j) stays valid for a virtual displacement from a equilibrium state to an arbitrary neighbouring state, 

dS" = p v dT» v - ^dN" + 2 - (A8) 

Here, the last term 2 denotes the contribution from the second order deviation from equilibrium. This postulate 
enables us to determine the form of the entropy flux in a near equilibrium system perturbatively. By addition of (|A7|) 
and (|A8p the the entropy four flux in the IS theory is obtained, 

= + T^p v - ^N" - Qf, (A9) 

where is an undetermined second order term in the deviations — T^ u , — Nq . 
We introduce, again, the energy-momentum tensor and particle flux as follows; 

T>* v = (s + P + n) - g" v (P + n) + vr^, (A10) 
AT" = nu" + iA (All) 

It should be noted that, different from the extended irreversible thermodynamics, the thermodynamic variables e, P 
and n satisfy the usual thermodynamic relations by construction. By using Eqs. (|A10p and (jAllj) . we finally obtain 
Eq. (jAU). 



APPENDIX B: INTERNAL- VARIABLE THEORY 

The internal- variable theories (IVT) is another approach to derive the generalized hydrodynamics (IH HH ■ Similarly 
to the extended irreversible thermodynamics, the IVT includes additional variables except for the usual thermody- 
namic variables. In this appendix, we use the idea of the IVT and rederive the IS theory. 

First, we assume the following modified first law, 

TdS = dE+ pdV - fidN - TdQ, (Bl) 

where Q is an additional variable. As we will see later, this definition of the modified first law is different from that 
in the extended irreversible thermodynamics. Then, the entropy production is given by 

°% ~ d » («"") = + ^ - ""V - a^. (B2) 

To obtain the same result as the IS theory, we assume Q as follows, 

Q = (AIT 2 + p 2 n^ + fou^) . (B3) 
Then, for the positivity of the r. h. s. of Eq. (|B2[) , we obtain the IS theory, 

dn nT d f0i\, a n, 

dr 2~ dr V ~T ] ~2~' 



d^u" + fa — + — — £ + Pi -flX ) = - ai U, (B4) 



P^ (d a up - p. 



dlT a /3 7T Q/3 T d f P 2 \ 7T Q/ 3 



dr 2 dr \ T J ' 2 

P" , (Vl+ l^ + ^w + l^^) ----- 



-A -FflX) = <*2lf", 1 Br,) 

The main difference between the extended irreversible thermodynamics and the IVT is the definition of the first 
law. By substituting Eq. (|B3j) into Eq. (|Blj) . we obtain 

TdS = dE + PdV - ndN - frndfi - AjTr^dTf^ - /?3^ M - Tn 2 d f - Tn^d (^) - Tv^d 

(B7) 
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One can easily see that the last three term on the r. h. s. do not exist in the extended irreversible thermodynamics. 
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